Spontaneous breaking of spatial and spin symmetry in spinor condensates 
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Parametric amplification of quantum fluctuations constitutes a fundamental mechanism for spon- 
taneous symmetry breaking. In our experiments, a spinor condensate acts as a parametric amplifier 
of spin modes, resulting in a twofold spontaneous breaking of spatial and spin symmetry in the 
amplified clouds. Our experiments permit a precise analysis of the amplification in specific spatial 
Bessel-like modes, allowing for the detailed understanding of the double symmetry breaking. On res- 
onances that create vortex-antivortex superpositions, we show that the cylindrical spatial symmetry 
is spontaneously broken, but phase squeezing prevents spin-symmetry breaking. If, however, non- 
degenerate spin modes contribute to the amplification, quantum interferences lead to spin-dependent 
density profiles and hence spontaneously-formed patterns in the longitudinal magnetization. 
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Spontaneous symmetry breaking is fundamental in dis- 
parate scenarios in physics ranging from cosmology [l| 
and particle physics [2[ to liquid crystals [Sl and su- 
perfluid Helium [4|. Symmetry breaking is also crucial 
in Bose-Einstein condensates (BECs), where the U{\) 
symmetry is spontaneously broken. Even more interest- 
ing, quantum gases provide unprecedented possibilities 
for the investigation of non-equilibrium dynamics, and 
in particular for the detailed analysis of dynamical sym- 
metry breaking, including the formation of topological 
defects via the Kibble-Zurek mechanism [l|, y, l5[- 

Spinor gases, formed by atoms with non-zero spin, 
present an exceptionally rich physics due to the non- 
trivial interplay between internal and external degrees 
of freedom [6|-|9[. A remarkable scenario is provided by 
BECs initially prepared in an unstable rap = (|0)) Zee- 
man state. Spin changing collisions lead to the correlated 
creation of atom pairs in mp = ±1 (|±1)) in a process 
resembling optical down conversion in nonlinear crystals. 
As a result, parametric amplification of matter waves in 
spinor BECs has been recently observed [IC 
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Parametric amplification of quantum fluctuations con- 
stitutes a fundamental mechanism for symmetry break- 
ing, and hence spinor BECs constitute an exciting sys- 
tem to investigate spontaneously broken symmetries. In 
a recent seminal experiment [5[ the sudden quench of 
a spinor BEC from a polar into a ferromagnetic phase 
was followed by the formation of ferromagnetic domains 
and topological defects in the transverse magnetization, 
whereas the longitudinal magnetization remained negli- 
gible. Ref. [5| provided a major insight in the formation 
of topological defects, but the nature of the symmetry 
breaking mechanism remained largely unexplored. 

This article explores the nature of the symmetry break- 
ing in a ^''Rb F = 2 BEC initially prepared in the |0) 
state. In our experimental system the |±1) atoms are 
amplified on well resolved magnetic field resonances [12| 
with characteristic Bessel-like spatial modes, since the 
energy splitting of these modes is on the same scale as 



the spin changing interactions. Our experiments hence 
permit a precise analysis of the amplification of quantum 
spin fluctuations, allowing for a detailed characterization 
of the symmetry breaking mechanism. Interestingly, a 
twofold spontaneous breaking of spatial and spin sym- 
metries may occur. On one hand, we show that quantum 
fluctuations of the relative phase between amplified de- 
generate vortex/antivortex states may break the cylin- 
drical symmetry imposed by the trap. On the other 
hand, contrary to the situation in Ref. [5|, the den- 
sity profiles in |±1) may differ from each other, lead- 
ing to spontaneously-formed longitudinal magnetization 
patterns only if various non-degenerate spin modes are 
significantly amplified. We show that this novel type of 
spin-symmetry breaking is linked to quantum interfer- 
ences occurring during the amplification. 

In the following we consider the initial BEC in |0) as 
a classical field yjn[){r), whereas the excitations in |±1) 
are represented by the operators 5ip±i{r) [l^l- If the 
population in |±1) remains small compared to that in |0), 
the dynamics may be described by the linear Hamiltonian 



H= d\J2 ^i'lHeffS^'. 
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where H^ff = -h'^\7^/2M + Veff{f) + q, M is the 
atomic mass, and VLeff{r) = Uino{r) characterizes the 
pair creation due to spin changing collisions. The effec- 
tive potential experienced by |±1) atoms is T4//(r) = 
F(f) + {Uo + Ui)no{r) — fi, with the harmonic trap F(r), 
and the chemical potential fi of the initial BEC. The in- 
teraction strengths are given by Uq and Ui [y, [ll[ . The 
magnetic field B influences the dynamics only due to the 
quadratic Zeeman energy, q cc B^ < 0. 

In our experiments the initial BECs are produced in 
a cylindrically symmetric trap well within the Thomas- 
Fermi regime. Hence V{r) + Uono{r) — fj, ~ holds within 
the BEC region r < rtf, where rtf is the Thomas- Fermi 
radius. Therefore Vef f is flat, apart from the small re- 
pulsion [/ino(r) which produces a bump of ~ /i x 30 Hz 



in the center. Beyond the BEC surface, Veff{r) = V{r) 
equals the external potential and rises sharply. Hence 
the radial potential may well be approximated by a 2D 
circular box potential with radius ~ r^f. In our experi- 
ments the energy splitting between transverse modes is 
larger than the typical energy of spin-changing collisions, 
and hence the transverse modes are well defined. On the 
contrary, the energy spacing of the accessible axial modes 
is much smaller than the spin-changing interaction, and 
hence they are not resolved. Our detection integrates 
along the axial direction, and hence a 2D description 
based on the transversal modes is sufficient. The single- 
atom eigenfunctions of the 2D circular box are 
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with cigcnencrgics e„; — fi^0^^j2rar1^. Here, Ji are Bessel 
functions of the first kind and /?„; is the nth zero of J; . 
Each mode is identified by two quantum numbers, n and 
Z, for the radial excitation and the angular momentum. 

A deeper understanding is gained by expanding ([T]) in 
these eigenfunctions: ^iimir) = Y.ni'^ni{r)an,i.m- The 
Hamiltonian then splits into H = J^r, 



\i\Hn\i\, where 



Hn\l\ - {in\l\ + q) X! "'n,l,mO-n,l,m 
7n,l=±\l\ 

+ ^ («i,i,l^n,-/,-l + "'i.-Ll^l,l.-l + ll-C-) . (3) 

[l3| . The mode is stable if the eigenen- 



with fl = UiUq 
ergy ^ni = \/(eri; + q)^ - ^"^ of Hni is real, whereas an 
imaginary ^„i leads to an exponentially increasing pop- 
ulation of ipni- Resonances in the spin dynamics occur 
when Eni + q ^ 0. Since e„i depend on \l\, modes with 
the same n and opposite / and —I are degenerate. 

The experimental apparatus to observe these Bcsscl 
modes and to study the symmetry breaking was de- 



scribed previously [IJ, llj|. Briefly, a BEC of ^"^Rb 



with a thermal fraction below 20% is prepared in the 
F = 2,mp ~ state in an optical dipole trap. Two 
intersecting laser beams at a wavelength of 1064 nm pro- 
vide a cylindrical trapping potential with measured oscil- 
lation frequencies of 187, 183 and 67 Hz. Before starting 
the experiment a strong magnetic field gradient is used 
to expel residual atoms in |to 7^ 0) states from the trap. 
This reduces unwanted atoms in |±1) which may domi- 



nate the amplification IjJ. Subsequently, the magnetic 
field is set to a value between and 2.5 G, and the spin 
components evolve for an adjustable hold time. 

After this spin evolution time, the trap is switched off 
to allow for ballistic expansion. During the first 2 ms, all 
components expand together, such that the |±1) clouds 
experience the repulsion from the expanding |0) BEC. 
Once the densities are low enough, the spin components 
are separated by a magnetic field gradient applied for 
3.5 ms. After a final expansion time of 1.5 ms, the 2D 
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FIG. 1. (a) Fraction of atoms transferred into the |±1) state 
within 18.5 ms as a function of the applied magnetic field. 
Each data point is an average over 30 realizations. The verti- 
cal gray lines indicate the resonance positions obtained from 
a 2D circular box model, and the labels indicate the corre- 
sponding Bessel inodes. (b) Absorption image of a [0) BEC 
and the |±1) clouds recorded at 1.29 G. 
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FIG. 2. Experimental and theoretical density distributions 
on the resonance positions after time-of- flight expansion, (a) 
Averaged experimental density profiles, (b) Calculated pure 
Bessel distributions corresponding to the experimental situ- 
ation, (c) Individual experimental density profiles, (d) Cal- 
culated superpositions of Bessel distributions (see text). The 
jo) BEC was omitted in (a) and (c) for clarity. 



column densities are imaged along the weak trap axis 
as shown in Fig. [T] (b) . Note that the expansion of the 
|±1) clouds is governed by two contributions. The ki- 
netic energy distribution in the trap and the repulsion 
from the |0) BEC during the initial expansion. Since 
the kinetic energy associated to the typical momentum 
~ l/^tf is negligible compared to the mean-field repul- 
sion, the clouds expand self-similarly [15|, ll6| . Therefore, 
the image obtained after time-of-flight expansion repro- 
duces the density distribution in the trap. 

Figure [T] shows the fraction of atoms in |±1) as a 
function of the magnetic field. One observes six reso- 
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FIG. 3. Relative position of the |±1) clouds prepared at 
0.64 G after variable hold time at high magnetic field. A 
fit with a damped oscillation yields a frequency of 40 ± 3 Hz, 
which closely corresponds to the energy difference between 
the resonances corresponding to the (1, 0) and (1, ±1) modes. 



nances, which arise when the quadratic Zeeman energy 
matches the eigenenergy of a mode (n, I). These positions 
are indicated by gray hnes, obtained from e„^i by fitting 
rtf = 3.9 /xm, in agreement with a calculated 3.7 /im. 
This is reflected by the density distributions of the |±1) 
clouds. Figure [2] (a) shows the averaged density pro- 
files of all realizations recorded at the six resonance po- 
sitions, which match the theoretical Bessel distributions 
of Fig. [2] (b) strikingly well. As expected the averaged 
profiles share the symmetries of the Hamiltonian ([1]), i.e. 
a spatial cylindrical symmetry and a spin symmetry un- 
der the interchange |1) <->■ | — 1). Only the resonance for 
(n, I) = (1, ±1), shows an orientation (see below). 

However, the individual density profiles shown in 
Fig. [21(c) are clearly different. Since the (1,0), (2,0) and 
(3, 0) modes have I = 0, the clouds remain cylindrically 
symmetric on the corresponding resonances even in in- 
dividual realizations. The remaining resonances however 
correspond to |/| = 1 modes, which are two-fold degen- 
erate, allowing for a clockwise rotation (vortex) or an 
anti-clockwise one (anti- vortex). Since the J±l) clouds 
result from amplified quantum fluctuations |ll| . the rel- 
ative strength of vortex/antivortex modes is chosen ran- 
domly for each realization. If just one of the modes is 
populated, the result is a cylindrically symmetric density 
distribution. However, if superpositions of the two modes 
are populated, non cylindrically-symmctric standing ra- 
dial waves emerge with a clear orientation (Fig. [2] (c) 
(2,±1) and (3,±1)). The angle of the orientation de- 
pends on the phase difference between the vortex and 
anti-vortex modes, which is randomly chosen by quan- 
tum fluctuations in individual realizations, constituting 
a clear example of spontaneous symmetry breaking. 

This argument does not hold for density profiles 
recorded on the (1,±1) resonance. Since it is very close 
to (1, 0) mode the clouds are produced in superpositions 
of (1,0) and (1,±1). As an example, a 50/50 superposi- 
tion is shown in Fig. [2] (c) and (d). Such superpositions 
are strongly off-centered and thus small magnetic field 
gradients are enough to break the cylindrical symmetry. 
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FIG. 4. Orientation of individual density distributions 
recorded in the vicinity of the (2, 1) mode, (a), (b) Distri- 
bution of angles for the |1) and | — 1) states, (c) Distribution 
of the difference between the angles of orientation. 



To test this interpretation, |±1) clouds were produced 
on this resonance and subsequently the magnetic field 
was set to a high value not supporting further spin dy- 
namics. After a variable hold time, the relative center of 
mass positions of the |±1) clouds was recorded. Figure[3] 
shows that the two clouds perform a damped oscillation 
at 40 ± 3 Hz, which matches the energy difference be- 
tween the / = and |/| = 1 resonances of 37 ± 5 Hz. This 
confirms the amplification of a coherent superposition in- 
cluding the (1,0) mode. 

Symmetry breaking is therefore best studied on the 
resonance corresponding to the (2,±1) mode. To eval- 
uate the symmetry breaking process statistically, 2000 
images were taken at four magnetic field values between 
1.58 G and 1.84 G for an evolution time of 21.5 ms. For 
each picture, the angle of orientation of the |±1) clouds 
was determined by two independent methods. First, a 
superposition of the two counter-rotating Bessel modes 
was fitted to the density distribution. Second, the eigen- 
vectors of the quadrupole moment tensor of the density 
distribution were determined. The difference of the two 
angle measurements has a standard deviation of 15° and 
78% of the measurements differing less than 40° were in- 
cluded. The remaining 22% were excluded from the anal- 
ysis, since the total atom number was too small or the 
image was cylindrically symmetric or experimental noise 
corrupted one of the angle measurements. We checked 
that the experimental result remained stable even for a 
large variation of the exclusion conditions. 

The results are presented in Fig. 2] For each magnetic 
field, a histogram shows the distribution of angles for 
|1) and | — 1) as well as the difference between these an- 
gles. All angles of the |±1) clouds occur with roughly the 
same probability, showing that the cylindrical symmetry 
is indeed spontaneously broken. This confirms that nei- 
ther the spurious production of |±1) atoms [ll| nor slight 



asymmetries of the potential influence the population of 
higher modes and thus do not break the symmetry [l7|. 

Interestingly, the relative angle (Fig. |4] (c)) shows an 
intriguing behavior. At 1.78 G, the two angles are locked 
to each other within the angle measurement uncertainty. 
Towards lower or higher magnetic fields, the locking re- 
laxes and differing angles become more probable. Such 
realizations with different orientations of the density dis- 
tributions in |1) and | — 1) break the spin symmetry of 
the density profiles. As a result, a spatial pattern of the 
longitudinal magnetization, iSz = X]m"^'^'/^m(^)'^V'm('^) 
emerges, contrary to Ref. [5| where Sz was negligible. 

This novel spin-symmetry breaking cannot be ex- 
plained only by the amplification of quantum fluctua- 
tions. The physical mechanism may be understood from 
Eq. ([3]), which for unstable modes may be re- written as 



Hn,\l\ — 
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with cos20„.; ~ {^n.\i\ + l)l^- The operators X^^^l\^^ 

can be identified as the quadratures ([-^„ ; ^'"^1 z nJ ~ 
i) corresponding to the amplification dynamics. The 



(2) 



quadrature X)^ [ ^ is stretched (~ e'^"!'! '*), whereas Xn^ 
is squeezed (~ e^l^"'i'il*). Hence the amplification of the 
(n, /)-mode is characterized by the mixing angle 0„.;. 

After a given time t we may evaluate the density- 
density correlation functions G,„.,„'(r) = {n„i{r)nmi {r)) 



G„,„(f) = 2Vi^a,n'(^), 



n.fl' 



(6) 
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J^n,«'(^,t)[l+COs2(0fl-fts')], (7) 



where only unstable modes n = (n, /) are summed, and 

p2(|5s|+|Cs,|)t 



Fs,s'(r) 



16 sin^ 26*^ sin^ 26*^, 



\Mm^n'{rr. (8) 



Equation ([7]) shows that spin symmetry breaking is 
caused by quantum interference between unstable modes 
with different mixing angles On- This is the case at mag- 
netic fields slightly above or below 1.78 G, for which the 
(2,0) or (3,0) modes interfere with the strong (2,±1) 
modes. Note, however, that degenerate modes (e„.i = 
en,-i) possess equal 0n,i-, resulting in the observed locked 
angles of the |±1) clouds exactly at 1.78 G. 

Alternatively, the angle locking may be explained by 
phase squeezing. At 1.78 G only the (2, ±1) modes 
contribute, and hence only ^2, 111 niust be considered. 



At the resonance (e2.|i| + (? = 0) the evolution op- 
erator e~'^^'iii*/'' corresponds to a two-mode squeez- 
ing operator, which produces two output waves with a 
perfectly anti-squeezed phase difference but a dynami- 
cally squeezed phase sum [18| when applied to the spin 
vacuum. For a sufficient evolution time, the phase 
sum approaches the squeezing angle [18|, and hence the 
phases 4>i=±i_m vn the linear superposition of (2,±1) ful- 
fill 4>i_i + 0_i__i = ^1,-1 -I- 4>-i,i- The orientation of 
the patterns is defined by the relative phase between the 
two counter-rotating modes in each \m) state, therefore 
the two angles are dynamically locked to each other since 
01,1 — 01,-1 = 0-1,1 ^ 0-1,-1- Again, the unlocking of 
the angles is due to the interference with neighboring 
resonances, which correspond to other Bessel modes. 

In summary, the amplification of |±1) states is char- 
acterized by a non-trivial double spontaneous symme- 
try breaking. On one hand, the cylindrical spatial sym- 
metry is broken by quantum fluctuations of the relative 
phase between vortex/antivortex modes, although phase 
squeezing may prevent the breaking of the spin symme- 
try. On the other hand, if various non-degenerate spin 
modes are amplified, spontaneously-formed longitudinal 
magnetization patterns appear due to quantum interfer- 
ences between the different amplified quadratures. Our 
results show that spinor gases constitute an exceptionally 
controllable system for the detailed analysis of symme- 
try breaking and in particular of its close connection to 
multimode squeezing during parametric amplification. 
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